Weakly interacting Bose gas in the one-dimensional limit 
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We prepare a chemically and thermally one-dimensional (Id) quantum degenerate Bose gas in a 
single microtrap. We introduce a new interferometric method to distinguish the quasicondensate 
fraction of the gas from the thermal cloud at finite temperature. We reach temperatures down to 
kT ~ 0.5hui± (transverse oscillator eigenfrequency u±) when collisional thermalization slows down 
as expected in Id. At the lowest temperatures the transverse momentum distribution exhibits a 
residual dependence on the line density nid, characteristic for Id systems. For very low densities 
the approach to the transverse single particle ground state is linear in nid- 
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Low dimensional systems are more fragile and prone 
to quantum and thermal fluctuations than their three- 
dimensional counterparts, so that the physics is drasti- 
cally changed [l[ . This is relevant in many different areas, 
but real systems are embedded in three dimensions. Cold 
atomic gases are well suited to synthesize low dimensional 
systems of interacting particles and to explore the influ- 
ence of the frozen dimensions. The dimensionality can 
be controlled by the external trapping potential, quan- 
tum degeneracy can be reached by established cooling 
techniques, and inter-particle interactions can be tuned 
by controlling the density or a Feshbach resonance Q • 

A number of experiments have been performed with 
two- (2d) and one-dimensional (Id) cold atom systems, 
observing phenomena like the Kosterlitz-Thoulcss tran- 
sition in 2d 0d^| or the 'fermionization' of bosons in the 
Id Tonks- Girardeau gas [1,0]. To date, deeply Id sys- 
tems have been studied mostly in optical lattices where 
many Id systems (typically ~ 100) are prepared at once 
@, 0]. These systems are then characterized by aver- 
ages over all 'tubes'. In contrast, creating only one or 
two Id clouds allows direct measurements of density Q 
and phase H, Til fluctuations, facilitates studies of dy- 
namical effects [l2| and allows to probe the full distri- 
bution function of quantum variables (l34l5j . A recent 



experiment has shown that the exact Yang- Yang theory 
predicts a density distribution of a single Id gas that 
matches the data better than ideal gas and mean field 
models (l6l | . Earlier experiments have been performed in 
elongated macroscopic magnetic traps, where the systems 
exhibited some characteristic Id properties while still in 
the ld-3d crossover (chemical potential fj, and tempera- 
ture kT > huj±). Most notably, phase fluctuating 'qua- 
sicondensates' were observed [17H19(. 

In this work we experimentally study the low tem- 



perature behavior of quantum degenerate Bose gases in 
the one-dimensional limit {n,kT < hujj_). We study the 
transverse expansion of an individual degenerate Id Bose 
gas, enabling us to explore the decreasing influence of 3d 
properties down to very low temperatures and chemi- 
cal potentials. When kT is lowered, the population of 
excited transverse modes is gradually reduced and the 
quasicondensate fraction grows. We directly investigate 
this process by observing the interference of two clouds 
after transverse expansion. The interference allows us 
to distinguish quasicondensed and thermal fractions of 
the conventional techniques fail in Id. We derive 

the temperature of the gas from the width ratio of these 
fractions and reach kT w 0.5?iu;_l as the lowest temper- 
ature in our experiment. Reaching such low T allows us 
to observe the exponential suppression of thermalization 
rates when the dimensionality is reduced from three to 
one. In the low temperature regime, we show that even 
for (i < hujj_, inter-particle interactions lead to an ad- 
mixture of 3d character at any finite density. Only in 
the limit of vanishing density, the ground state width 
expected for the noninteracting gas is recovered. 

In our experiment we prepare Id Bose gases of 87 Rb 
in a microscopic magnetic trap formed at a distance of 
d ~ 10 — 30 /im from the surface of an atom chip 2^, 21 1 
as described in detail in [22] |. At small atom-surface dis- 
tances, small longitudinal variations of the trapping po- 
tential arising from slightly non uniform currents in trap- 
ping wires modulate the atomic density of the quantum 
degenerate gas. It is essential that these modulations 
are smaller than the chemical potential of the cloud and 
allow us to study linear density nid dependent effects 
[23I ] . We are able to explore the entire density range from 
the detection threshold of ~ 3^m _1 , correspond- 
ing to an interaction energy per particle of the order of 
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FIG. 1: (color online) Interference of two Id quasiconden- 
sates at relatively low (top left) and high (bottom left) tem- 
peratures. The absorption images were taken along the long 
direction (z) of the clouds after they had expanded from a 
double well trap in potential free time-of-fhght. The absolute 
value of the Fourier transform (in arbitrary units) along the 
direction (x) of modulated density (fringes) is also depicted. 
The central panels show the density profiles at the interference 
fringe period (dashed blue line) and at zero spatial frequency 
(solid red line) and Gaussian fits (bars indicate the respective 
standard deviations a). The panel on the right illustrates the 
increase in ratio between the widths of the cloud and the part 
of it that interferes as a function of final trap depth after the 
evaporative cooling process. 



0.03 hu± only (iv± = 2tt x 4kHz), to well into the ld-3d 
crossover regime /x > hui± {nid w 250yum~ 1 ). The length 
of the continuous gas reaches L sa 1 mm, corresponding 
to an aspect ratio of L/a± > 5000 (a± = 170 nm is the 
transverse ground state size) , the Lieb-Liniger parameter 
ranges between 7 = 0.0014 — 0.12. 

At finite temperature, an important question is which 
fraction of the gas is (quasi)condensed. The time-of-fught 
(TOF) technique of observation (in free expansion) of bi- 
modal density profiles with a parabolic core and a sur- 
rounding thermal cloud has been very successful in stud- 
ies of 3d BECs. In the case of Id gases, however, this 
approach cannot be easily applied since both quasicon- 
densate and thermal fraction expand transversely in a 
Gaussian fashion. Their widths after expansion hardly 
differ at sufficiently low T. 

Our solution to this problem is to exploit the coher- 
ence property of a quasicondensate as revealed in an in- 
terference experiment. We trap a gas at sufficiently large 
distance from the surface (d ~ 30 (im) to obtain com- 
pletely smooth Id clouds. We then split the gas using a 
radio frequency induced transformation of the trap into 
a double well as described in [2414261]. These two clouds 
are then released from the confining potential by sud- 
denly switching off all trapping fields. After a free TOF 
expansion of 16 ms, the clouds overlap and form a high 
contrast interference pattern wherever the single particle 
ground state is macroscopically occupied. 

Fig. Q] shows two example images obtained in the way 
described above. One represents the low T regime where 
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FIG. 2: (color online) Collision rate devided by density as 
a function of temperature across the ld-3d crossover. The 
lowest achieved temperatures (kT < 0.5huJ±) appear to be 
limited by slow rethermalization due to inhibited elastic col- 
lisions. The increase in collision rate as temperature is in- 
creased does not follow the square root law (solid line shows 
best fit) that is valid in 3d, but is compatible with an expo- 
nential suppression as predicted in Id (dashed line shows best 
fit for data within kT < 2hio±) [zj]]. 



essentially the entire cloud interferes, i.e. the sample is 
close to a pure quasicondensate. The other image clearly 
displays a core that interferes with high contrast, over- 
lapped by a larger structureless thermal cloud. For quan- 
titative analysis we perform a one-dimensional Fourier 
transform (FT) along each pixel row of each image. We 
determine the width of the coherent quasicondensate 
fraction by fitting a Gaussian to the first harmonic at fen, 
the wave vector corresponding to the period of the fringe 
pattern. For comparison we fit a Gaussian to the over- 
all density profile. Note that the conventional method of 
fitting bimodal profiles fails already at kT rj 3Huj±. 

The total transverse- momentum distribution of the ex- 
panded cloud is wider than its coherent part since the 
excited states of the radial motion contribute to the for- 
mer, whereas the latter is determined by the properties 
of the ground state of the radial motion. The experi- 
mental parameter controlling the temperature of the gas 
is the final frequency i/ r [ of the radio frequency (rf) field 
(measured from the trap bottom) used to reduce the trap 
depth in the applied forced evaporation scheme. The ra- 
tio between thermal and coherent parts of the cloud as 
a function of the final depth of the trap hvrf during the 
cooling process is plotted in Fig. [1] (right). The coher- 
ent cloud's width indeed does not depend on z/ r f, small 
remaining variations (see also Fig. [3]) can be attributed 
to variations in atom density as discussed below in more 
detail. At the lowest temperatures, the trap depth is set 
to very close to the chemical potential of the gas and 
partially depletes the quasicondensate. At high temper- 
ature, the quasicondensate fraction decreases as the crit- 
ical degeneracy temperature is approached. In contrast, 
the width of the thermal cloud increases with his r f, as 
expected. 
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We assign a temperature T to the gas according to 
the transverse (y) width a y expected from an ideal Bosc- 
Einstcin gas. The temperature dependent width upon 
expansion was obtained by calculating the expectation 
value (y 2 ) using the grand canonical partition function 
of a 3d non-interacting Bose gas. Reflecting the experi- 
mental implementation, this gas is strongly harmonically 
confined in the two transverse dimensions. The chemical 
potential \i was fixed by the particle density at each T. 
We have verified that our experimental parameters are 
well in the dilute gas limit where interactions do not lead 
to detectable deviations from the ideal gas behavior [2|| . 
The lowest temperature we can reliably measure by this 
procedure is kT w 0.5?iu;^. 

In evaporative cooling the ratio of hv r { to the final 
temperature kT (truncation parameter 77 — hv^/kT) is 
a measure of the efficiency of the cooling. For kT > huj±_ 
we find 77 ~ 5. As we approach colder temperatures 
kT ~ ?ia;_L, 77 drops significantly and reaches rj rj 1 for the 
lowest measured T. In our experiments the atoms were 
always released after a constant hold time t = 100 ms in 
the final trap. This way 77 can be used to quantify elas- 
tic collision rates j c leading to thcrmalization of the gas 
[1^. Fig. [5] displays 7 C as a function of T. In 3d 7 C is 
expected to be given by the constant elastic collisional 
cross section, the density n, and the mean particle veloc- 
ity, i.e. 7c/^ oc VT. The data clearly deviate from this 
scaling behavior. The faster drop of 7 C for low T can be 
related to reduced thermalizing binary collisions in the Id 
regime. To illustrate the possibility of this effect being 
relevant in our case, we have fitted two model functions 
to our data: 1.) A square root behavior as expected in 
3d; 2.) an exponential drop of thermalizing binary colli- 
sion rates according to 7 c /tt. oc cxp (— 2huj±/kT) [27| as 
expected in the low temperature regime (kT < hujj_ , the 
fit is for kT < 2hu)±). Three-body collisions start to be 
the dominant mechanism for kT < 0.37kdj_ [27, 3^, 3l| . 
We find that the strong reductions of thermalizing elastic 
binary collisions is consistent with our data. The associ- 
ated reduction of thcrmalization rates 32j is likely to be 
the limiting factor in cooling further into the Id regime. 

We now investigate the low temperature limit ksT < 
Hluj_ for a single Id quasicondensate, in order to uncover 
the density dependent expanded cloud widths that can- 
not be seen at higher temperatures 33j . Here we exploit 
the density modulations that occur for small trap-surface 
distances (d = 10 71m). The inset of Figure |3] shows a 
typical density profile of an expanding cloud after 6 ms 
TOF in a two-dimensional projection onto the xz plane, 
so that the longitudinal (2) and one transverse (y) direc- 
tions are resolved. The density profile is strongly modu- 
lated along z. The variations on longer length scales arise 
from a correspondingly inhomogeneous (and stable) po- 
tential. As the gas does not notably expand along the 
weakly confining z direction during TOF, these features 
of the density profile are unaffected by TOF. We can con- 
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FIG. 3: (color online) Transverse width (Gaussian standard 
deviation, normalized to the single-particle ground state ex- 
panded size) of the expanded cloud after time-of-flight as a 
function of longitudinal line density of atoms throughout the 
quasi one-dimensional regime of /i <C hcu± and into the ld-3d 
crossover. The points represent our data, measured in two 
different traps (blue squares and red circles, example data 
see upper inset). The line shows the theory without any 
adjustable parameters (Eq. [T] HU). The lower insets illus- 
trate that a Gaussian distribution models the data very well 
in the Id regime (left, nid = 100/im -1 ), while systematic 
deviations become discernable in the ld-3d crossover regime 
(right, nid = 250/^m -1 ). Data, best Gaussian fits and resid- 
uals (magnified by a factor of 3) are shown. 



firm this by comparison to in situ images taken under 
the same experimental conditions. The known quasicon- 
densate nature of the Id quantum degenerate Bose gas 
causes random phase fluctuations in the trapped system 
that are converted to short range density fluctuations af- 
ter TOF 0, 
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18j . This randomness allows us to 
eliminate this type of fluctuations by averaging over a few 
TOF images to ensure that the remaining modulation is 
due to actual density modulation in the trapped gas. In 
the transverse y direction, on the other hand, the cloud's 
initial extension is comparable to the harmonic oscillator 
ground state size of ~ 170 nm and cannot be resolved by 
the in situ imaging. After TOF, the size of the cloud 
along y is then a direct measure of the energy stored in 
the trapped gas, since no longitudinal expansion occurs 
on the experimentally relevant timescales. 

In the common case of a 3d harmonic trap when /i 3> 
Hluj_, the interaction energy dominates the kinetic energy 
and a parabolic density profile emerges in the Thomas- 
Fermi approximation. In contrast, the zero point energy 
hujj_/2 dominates in the Id case when fi -C hu>x, so that 
the shape of the distribution becomes Gaussian. By ap- 
plying a local density approximation to averaged sets of 
images of the type shown in the inset of Fig. [3J we ob- 
tain transverse TOF profiles of the cloud distribution as 
a function of n^- We compare these distributions to best 
fits to Gaussian distributions (Fig. [3] insets) and confirm 



4 



that as long as riid < 100 ^m^ 1 , the Gaussian fits rep- 
resent the data much more closely than the parabolic 
ones that would be expected in 3d. For larger densities 
nid ^ 100 /^m _1 , in the ld-3d crossover, where \i > htu±, 
systematic deviations from the fitted model become dis- 
cernable. Hence, Gaussian expansion is characteristic of 
the quasi Id regime in which the shape of the momentum 
distribution of the single particle transverse ground state 
is directly observed. 

Fig. [3] shows the width of the fitted Gaussian distri- 
butions as a function of nid- Even down to the low- 
est linear densities, repulsive inter-particle interactions 
(scattering length a scat = 5.3 nm) broaden the distri- 
bution. Note that the scattering is still of 3d nature 
(flscat "C aj_ = \J h / mjjj ±) , even in the regime where bi- 
nary collisions cease to contribute to thermalization when 
kT <C fi,w_i_.We compare our measured widths cr(nid) to 
the theoretical predictions of the correction to the width 
do of the single-particle ground state case [13, HBJ 

= (l + 4a scat n ld )* (1) 

^0 

and find excellent agreement in the absence of any fit- 
ting parameters for /i < tku±. The data only starts to 
slightly deviate form this model when \i sa hcj± , i.e. when 
the expansion of the gas is no longer purely Gaussian 
in the ld-3d crossover regime. For large temperatures 
(T ~ Tiwj_ or larger) the measured width saturates and 
stays at a constant value [33j given by the temperature, 
significantly above the curve shown in Fig. [3] 

In conclusion, we have prepared quantum degenerate 
bosonic gases in the one dimensional limit. We use an 
interferometric technique to differentiate between qua- 
sicondensed and thermal part of the gas at finite tem- 
peratures. We find that the width of only the thermal 
cloud grows with temperature while the coherent part is 
temperature independent. Determining the temperature 
from the thermal cloud width and comparing it to the 
trap depth indicates a loss of cooling efficiency at very 
low temperatures kT < huj±, which can be attributed 
to suppressed thermalization in the Id regime. In the 
low temperature regime (kT ps 0.57ia;_i_), the transverse 
momentum distribution exhibits Gaussian shape. The 
width of this distribution, however, equals that of the 
single particle ground state only in the limit of vanish- 
ing line density when 7 — > 00. This implies that a true 
Id gas is always a strongly correlated Tonks- Girardeau 
gas 0, 0]. Our measured widths follow the theoretical 
expectation for repulsive interactions. 
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